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ABSTRACT 

We present a Chandra X-ray Observatory ACIS-S variability, spectroscopy, and imaging study of 
the peculiar binary containing the millisecond pulsar J1023+0038. The X-ray emission from the 
system exhibits highly significant (12.5a) large-amplitude (factor of 2-3) orbital variability over the 
five consecutive orbits covered by the observation, with a pronounced decline in the flux at all energies 
at superior conjunction. This can be naturally explained by a partial geometric occultation by the 
secondary star of an X-ray-emitting intrabinary shock, produced by the interaction of outflows from 
the two stars. The depth and duration of the eclipse imply that the intrabinary shock is localized near 
or at the surface of the companion star and close to the inner Lagrangian point. The energetics of 
the shock favor a magnetically dominated pulsar wind that is focused into the orbital plane, requiring 
close alignment of the pulsar spin and orbital angular momentum axes. The X-ray spectrum consists 
of a dominant non-thermal component and at least one thermal component, likely originating from the 
heated pulsar polar caps, although a portion of this emission may be from an optically-thin "corona" . 
We find no evidence for extended emission due to a pulsar wind nebula or bow shock down to a limiting 
luminosity of Lx < 3.6 x 10 29 ergs s _1 (0.3-8 keV), <7 x 10~ 6 of the pulsar spin-down luminosity, 
for a distance of 1.3 kpc and an assumed power-law spectrum with photon index F = 1.5. 
Subject headings: pulsars: general — pulsars: individual (PSR J1023+0038) — stars: neutron - 
X-rays: stars 



1. INTRODUCTION 

The Galac tic source FIRST J102347.6+003841 was 
discovered bv lBond et al.l (|2002f ) who initially classified it 
as a magnetic cataclysmic variable . Subsequent studies 
(jThornstensen fc Armstrong! 120051 : iHomer et ail I2006T ) 
concluded that this object is more likely a neutron star 
low-mass X-ray binary. Its true nature was finally un- 
veiled with the Green Bank Telescope discovery of 1.69- 
ms radio pulsations (Archibald et al. 2009), establishing 
it as a compact binary containing a rotation-powered mil- 
lisecond pulsar (MSP). The pulsar, PSR J1023+0038, 
is in a circular 4.8-hour binary orbit with an optically 
identified non-degenerate low-m ass (~0.2 M ) compan- 
ion star (|Archibald et al.l l2009f ) . Perhaps more impor- 
tantly, it is the first radio MSP binary to have exhibited 
past evidence for an accretion disk. In 2001, double- 
peaked emission lines characteristic of accretion disks, ac- 
companied by short-timescale flickering and a blue spec- 
trum, were observed at op tical wavelengths ()Bond et al.l 
120021: iSzkodv et all 120031 ) . There is no e vidence that 
a dis k is currently present in the system (|Wang et al.l 
2009). A plausible interpretation of this behavior is a 
very recent transition from an accretion- to a rotation- 
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powered neutron star. If this is indeed the case, PSR 
J1023+0038 lends strong observational support for the 
long-suspected evolutionary connection between accret- 
ing MSPs in low-mass X-ray binaries ((Wiinands; I2010L 
and re ferences therein) a nd rotation-powered "recycled" 
MSPs (|Ah>ar et al.lll982D . 

In its present "quiescent" state, PSR J1023+0038 ex- 
hibits radio eclipses at superior conjunction, when the 
secondary star is between the pulsar and observer, as well 
as short, random and irregular eclipse s and dispersion 
meas ure variations at all orbital phases ((Archibald et al.l 
2009). Based on the orbital parameters of the binary 
and an assumed M = 1.4 M Q pulsar, for the implied 
orbital inclination of w46°, the line-of-sight between the 
pulsar and the Earth does not intersect t he Roche lobe of 
the c ompanion at any point in the orbit (|Archibald et all 
l2010f h Therefore, the eclipses must be caused by mate- 
rial driven off the surface of the secondary by the im- 
pinging pulsar wind. Such eclipses and dispersion mea- 
sure variations seen in si milar MSP binaries both in 
globular clusters (see, e.g.. iD'Amico et aT1l200lL for the 
case of PSR J1740 -5340 in NGC 6397 ) and in the field 
of the Galaxy fsee iHessels et all 120111 concerning PSR 
J22 15+51) can be attributed to the presence of matter 
flowing out from the irradiated companion, as is likely 
the case for J1023+0038. 

The J1023+0038 system has been previously ob served 
in X-rays with XMM-Newton (jHomer et al.ll2006f ) for 15 
ks, with all instruments in full imaging mode, and revis- 
ited in a target of opportunity observation for 34.5 ks, 
with t he EPIC pn instrumen t configured in fast timing 
mode (jArchibald et al.l l2010h . The latter observation, 
in addition to detecting variability correlated with or- 
bital phase, has revealed compelling evidence for X-ray 
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pulsations at the pulsar's rotational period. The X-ray 
spectrum is relatively hard with power-law spectral pho- 
ton index r « 1 — 1.3 and luminosity Lx « 9 x 10 31 
ergs s _1 (0.5-10 keV), plus a possible thermal compo- 
nent. The nature of the J1023+0038 system and the 
evidence for variability in the X-ray lightcurve favor in- 
teraction of the relativistic particle wind from the pul- 
sar with matter from the close stellar companion as 
the source of the non-ther mal radiation. The de tec- 
tion of X-ray pulsations by lArchibald et al.l (|2010f ) in 
the XMM-Newton EPIC pn data implies that a portion 
of the observed X-rays originates from the pulsar itself. 
Pulsed X-ray emission from typical MSPs (E « 10 33 ergs 
s" 1 ) is thermal in nature (iZavlirj [20061 : iBogdanov et al.l 
20061 IBoedanov fc Gr indlav 2009), while in the less com- 
mon energetic MSPs (E m 10 36 ergs s^ 1 ) the emission 
is characteristically non-thermal (see, e.g., Rutledge et 
al. 2004). Given that PSR J1023+0038 is moderately 
energetic, with E — 5 x 10 34 ergs s -1 , the true natu re of 
the pulsations is ambiguous (lArchibald et al.ll20Toh . 

Most of the rotational energy lost by a pulsar is car- 
ried away by win ds comprised of relativisti c particles and 
magnetic fields (jKennel fc CoronTtl Il984h . The winds 
shocked in the ambient medium produce pulsar wind 
nebulae (PWNe) observable from the radio through j- 
rays. To date, extended X-ray emission has only been 
observed from two MSPs: the canonica l "black widow" 
PSR B1957+20 (|Fruchter et al.l [1981 IStappers et all 
l200l and the nearby isolated PSR J2124-3358 (Hui 
& Becker 2006). As PSR J1023+0038 moves through 
the ambient interstellar medium, it may also generate 
sufficient ram pressure to confine the pulsar wind, re- 
sulting in the formation of a bow shock, accompanied 
by a synchrotron-emitting "cometary" tail trailing be- 
hind the pulsar. Given the strong possibility that PSR 
J1023+0038 is a recently activated radio MSP, it is im- 
portant to search for an associated wind nebula. The 
spatial structure of any X-ray tail along the proper mo- 
tion direction could, in principle, be employed as a "time- 
line" of the pulsar wind activity. This is possible because 
any accretion onto a radio pulsar would presumably ex- 
tinguish the pulsar wind, which in turn would result in 
interruptions in the injection of the energetic wind into 
the interstellar medium. 

Herein, we report on a Chandra X-ray Observatory 
Advanced CCD Imaging Spectrometer (ACIS) imaging 
spectroscopic observation of PSR J1023+0038. This in- 
vestigation offers further insight into the properties of 
this remarkable binary. The present work is outlined as 
follows. In §2, we summarize the observations, data re- 
duction and analysis procedure. In §3, we focus on the 
X-ray variability of PSR J1023+0038. In §4 we present 
the spectroscopic analysis, while in §5 we conduct an 
imaging analysis. In §6 we describe a simple geomet- 
ric model and in §7 the physics of the intrabinary shock 
emission. We offer conclusions in §8. 

2. OBSERVATION AND DATA REDUCTION 

The Chandra dataset was acquired on 2010 March 24 
(ObsID 11075) in a single, uninterrupted 83.1-ks expo- 
sure. The pulsar was placed at the nominal aim point 
of the back-illuminated ACIS-S3 CCD, configured in 
VFAINT telemetry mode. To minimize the effect of pho- 
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Fig. 1— Chandra ACIS-S lightcurve of PSR J1023+0038 in 
the 0.3-8 keV (Top panel) and 0.3-1.9 keV and 1.9-8 keV bands 
(Bottom panel). The five minima occur at </>(, fa 0.1 — 0.4 for 
each of the five full orbits covered by the observation. The orbital 
phase is defined based on the radio pulsar timing convention, in 
which superior conjunction (when the secondary star is between 
the pulsar and observer) occurs at </>), = 0.25. The background 
contributes only 0.1% to the total count rate. 

ton pileup (Davis 2001), the detector was used in a cus- 
tom subarray mode with 256 pixel rows, starting from 
CCD row 385. 

The data re-processing, reduction, and analysis were 
performed using CIAC0 4.2 and the corresponding cali- 
bration products (CALDB 4.3.1). Starting from the level 
1 data products of the ACIS-S observation, we first re- 
moved pixel randomization from the standard pipeline 
processing to aid in the search for extended emission near 
the point-source emission from the pulsar. To search for 
large-scale diffuse emission, we applied the background 
cleaning algorithm applicable to the VFAINT telemetry 
mode. However, this procedure tends to reject real source 
counts for relatively bright sources. Consequently, for the 
image subtraction, spectroscopic, and variability analy- 
ses discussed below, the data without the background 
cleaning applied were used. 

For the spectroscopic and variability analyses, we ex- 
tracted the emission from a circular region of radius 2" 
that encloses >90% of the total source energy at 1.5 keV. 
The pulsar net count rate in this region is 0.0394±0.0007 
counts s _1 (0.3-8 keV). To permit spectral fitting in 
XSPEC0 12.6.0q, the extracted source counts in the 0.3- 
8 keV range were grouped in energy bins so as to ensure 
at least 15 counts per bin. The background was taken 
from three source-free regions in the image around the 
pulsar. For the variability analysis, the photon arrival 
times were translated to the solar system barycenter us- 

8 Chandra Interactive Analysis of Observations, available at 

|http : //cxc .harvard ■ edu/ ciao/| 

a Available at http : //heasarc .nasa. gov/docs/xanadu/xspec/index .html 
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4. SPECTROSCOPY 
4.1. Total Spectrum 
As found by iHomer et all (|2006f) and lArchibald et al.l 



Fig. 2. — (Top panel) Chandra ACIS-S emission in the 0.3-8 keV 
range from PSR J1023+0038 folded at the 4.8-hour binary period. 
The dotted magenta line is the lightcurve smoothed using a four- 
point moving average. (Middle panel) Lightcurves in the 0.3—1.9 
keV (red) and 1.9—8 keV (blue) bands. (Bottom panel) Ratio of 
the lightcurves in the middle panel, with the 1.9—8 keV divided by 
the 0.3-1.9 keV count rate. For reference, the grey bands represent 
the approximate portions of the orbit where the pulsar undergoes 
radio eclipses at 0.35, 0.7, 1.6, and 3 GHz (from lightest to darkest 
grey, respectively). Two cycles are shown for clarity. 



ing the CIAO tool axbary assuming the DE405 JPL solar 
system ephe meris and the pulsar p osition derived from 
radio timing ([Archibald et al.l [20 09). We note that the 
3.2-s time resolution of the ACIS-S precludes the study 
of pulsations at the pulsar spin period. 

3. X-RAY ORBITAL VARIABILITY 



The study by Archib ald et al.l ([2010D uncovered large 
amplitude X-ray variability from J1023+0038. However, 
due to the short exposure, it was not possible to es- 
tablish whether the flux modulations were truly peri- 
odic. We ha ve used the radio timi ng ephemeris of PSR 
J1023+0038 ([Archibald et al.ll2009h to fold the barycen- 
tered Chandra source photons at the binary period. As 
the observation covers over five binary orbits, we can 
definitively confirm that the X-ray flux is modulated at 
the binary period (Figures 1 and 2). The variations are 
characterized by a factor of ~2-3 decline in the photon 
coun t rate at super ior conjunction ((f>i, w 0.25). A Kuiper 
test ([Paltanil [2004) on the folded, unbinned lightcurve, 
weighted to account for the non-uniform exposure across 
the orbit, indicates a 6 x 10~ 36 (12.5er) probability that 
photons being drawn from a constant distribution would 
exhibit this level of non-uniformity. Statistically, the in- 
dividual lightcurves from the five binary orbits are consis- 
tent with having the same shape. There is no indication 
of statistically significant spectral variability throughout 
the orbit (bottom panel of Figure 2). As discussed in 
§6, the X-ray modulations can be plausibly interpreted 
as being due to a partial geometric occultation of an in- 
trabinary shock by the secondary star. 



(2010) based on the set of XMM-Newton observations, 
the phase-integrated Chandra X-ray spectrum of PSR 
J 1023+0038 is well described by a pure absorbed power- 
law, while it is poorly fitted by a single thermal (black- 
body or neutron star hydroge n atmosphere) model. 
Moreover. lArchibald et all (|2010f ) found that a slight im- 
provement in the fit quality was obtained with the ad- 
dition of a thermal (neutron star atmosphere) compo- 
nent. Based on this, in the spectral analysis of the Chan- 
dra data we consider both a one-component absorbed 
power-law model and a composite absorbed power-law 
plus NS atm osphere. We use th e hydrogen atmosphere 
NSA model (|Zavlin et al.l [19961 over a blackbody be- 
cause an atmospheric layer is expected at the surface of 
a MSP given the standard formation scenario of MSPs, 
involving accretion of a substantial amount of matter. 
In addition, the thermal pulsations from the nearest 
known MS Ps require the presence of a light-element at- 
mosp here (Bogdan ov et al.l l2007t iBogdanov fc Grindlavl 
2009) . In principle, a portion of the observed X-ray emis- 
sion from the J1023+0038 system could also originate 
from a thermal plasma within or around the binary, pos- 
sibly from the active corona of the secondary star or the 
plasma responsible for the radio eclipses. To explore this 
possibility, we also conduct fits with a model consist- 
ing of a power-law plus MEKAL thermal plasma model 
in XSPEC. The MEKAL model considers an emission 
spectrum from hot diffuse gas including line emissions 
from several elem e nts g i ven a set of metal abundances 
([Mewe et al.lfl985l fl98l ILiedahl et al.1119951) . The best- 
fit parameters from the three different models are sum- 
marized in column 2 of Table 1 . 

Both the pure power-law and power-law plus thermal 
component yield statistically acceptable fits. However, 
an F-tcst indicates that the probability the composite 
model is a better fit purely by chance is only 1.2 x 10~ 6 , 
suggesting that the addition of the thermal component 
is warranted by the data. By analogy with most MSPs 
detected in X-ravs (iZavlinl [20061 : I Bogdanov etall 120061: 
IBogdanov fc Grindlavl 120091 ) . PSR J1023+0038 should 
have heated magnetic polar caps and hence a thermal 
component in the predominantly non-thermal spectrum. 
The evidence for re latively broad pulsatio ns at the pul- 
sar period found by lArchibald et al.l (|2010D in the XMM- 
Newton EPIC pn data are consistent with a portion of 
the observed X-rays originating from the neutron star 
polar caps. The implied thermal luminosity of 6 x 10 30 
ergs s _1 and E to Lx conversion efficie ncy (^10~ 4 
are comparable to those of most nearby (jZavlinl 1200 
IBogdanov fc Grindlav l 120091) and globular cluster MSPs 
(IBogdanov et al.ll2006l 120111) . 

As in lArchibald et al. nf2oioh . in most of the fits, the 
hydrogen column density along the line of sight is con- 
sistent with zero. This is not surprising, given that the 
column density along the line of sight and to the edge of 
the Galaxy is only Ah ~ 1 x 10 20 cm~ 2 . Since such a 
low Ah only affects the spectrum below ~0.3 keV, the 
fits of the Chandra spectrum are nearly insensitive to 
this parameter. 
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TABLE 1 

Summary of Spectral fits for PSR J1023+0038. 



Model 3 


Total 


06,1 

(0.0 - 0.5) 


4>b,2 

(0.5 - 1.0) 


Joint 

<t>l + 4>2 


Power-law 


N n (10' M cm" 2 ) 

r 

Fx (10 -13 ergs cm -2 s _1 ) c 
Xl/dof 


< 0.15 
1 1Q+ 003 

3.»U_ Q9 
0.99/163 


< 0.23 

1 2 n+0.05 

1.19/71 


< 0.48 

I 14+004 
1 ' 1 ^-0.03 

4 66+0- 15 
^• DD -0.12 

1.20/102 


< 0.17 

, 2O+0.05/-, , 4 +0.04 
2 77+°- 11 /4 66+ ' 15 
1.19/174 


Power-law 4- H atmosphere (NSA) b 


N H (10 20 cm- 2 ) 

r 

T cff (10 6 K) 

Reft (km) 

Thermal fraction d 

F x (10 -13 ergs cm" 2 s" 1 ) 

Xl/dof 


< 5.4 
1 00+ 05 
76+ ' 18 

06+ ' 04 

| ie + 0.15 
^• lo -0.12 

0.84/161 


1 09+ 07 
u - 44 -0.07 

16.3ll 8 6 

21 +°' 09 
u - zl -0.08 
o 71 +0.30 
' 1 -0.28 

0.85/69 


< 7.9 
1 oo+o- 10 

1 1+ 4 ' 4 
06+ 06 

4 07+O.21 
^• 9 '-0.21 

1.18/100 


7 3+ 2 - 4 

i-o»±S:S?7i.m±8:8? 

55+0- 09 
u -°°-0.06 

8.3l»;? 
0.17+g- g 7 / < 0.06 
3.49to ; 2 i/4.86lO;i 3 
1.06/172 


Power-law 4- MEKAL" 


Nh (10 2U cm" 2 ) 

r 

kT (keV) 

Thermal fraction d 

Fx (10 -13 ergs cm -2 s _1 ) c 

Xl/dol 


< 0.54 

, n7 +0.04 

>-■'•>< -0.03 

22+ ' 18 

n n oo+0.033 
U.UJO_ 032 

4 nq+ 09 
0.84/161 


< 0.72 

i 11+0.07 

J - li -0.05 

iq+ 02 
u - ly -0.03 

o.06±g;gg 
3.0310;^ 
0.86/69 


< 1.3 

1 - uo -0.05 

26+ 05 
u - ZD -0.05 

< 0.07 

A O1+0.15 

4 - Si -0.15 

1.16/100 


< 0.67 

1 1q +0. 06 /1 n ,+0.04 
1 - LO -0.06/ i - u,3 -0.04 

n 99 +0.02 

u ' zz -0.02 

0.17lO;0 7 /o.lllO;05 
2.97l0;0 9 /4.86l0; 43 
1.06/172 



a All uncertainties and upper limits quoted are la. 

b For the NSA model, a M = 1.4 Mq, i? = 10 km neutron star is assumed. The effective emission radius as measured at the 
stellar surface, R c ft, was calculated assuming a distance of 1.3 kpc (|Archibald et al.ll2009D . 
c Unabsorbed X-ray flux (0.3-8 keV) in units of 10 -13 ergs cm -2 s _1 . 
d Fraction of unabsorbed flux from the thermal component in the 0.3-8 keV band. 



(Figure 3). We fit the two spectra both separately and 
jointly. In the joint fits, we tied Nh and both the temper- 
ature and effective radius of the thermal component in all 
cases since the radiation from the pulsar is not expected 
to exhibit any dependence on orbital phase. The results 
of the orbital phase- resolved spectral fits are summarized 
in columns 3-5 of Table 1. 

Although statistically a pure power-law spectrum pro- 
vides a good description of the total phase-averaged X- 
ray spectrum, it yields only marginally acceptable fits 
of the individual spectra at orbital phases — 0.5 and 
0.5 — 1, with xt ~ 1-2 in both cases. The same is true 
for the joint fit of the two phase intervals. In the indi- 
vidual and joint fits, the inferred spectral photon indices 
are consistent with being the same. This is in agreement 
with the finding in §3 regarding the absence of orbital 
phase-dependent spectral variability although deeper ob- 
servations are necessary to confirm this. 

In general, the spectral fits tend to favor the non- 
thermal plus thermal model over the purely non-thermal 
model. It should be noted, however, that the derived 
temperatures and radii for the NSA model are consid- 
erably different between the total, individual and joint 
spectral fits. One possible explanation is the presence 
of a multi-temperature thermal spectrum, arising due to 
non-uniform heating across the face of the magnetic po- 
lar cap s. Such thermal emission is seen in several nearby 
MSPs (|Zavlidl2006t iBogdanov fc Grindlavl f2009h . The 
variation in the dominant non-thermal component as a 
function of orbital phase could result in the fits being 
more sensitive to the different temperatures such that at 



c For the MEKAL model, solar abundances are assumed. 
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Fig. 3.— Chandra ACIS-S X-ray spectra of PSR J1023+0038 for 
orbital phases cij> = 0.0 — 0.5 (black) and <j>b = 0.5 — 1.0 (red), fitted 
with a model consisting of an absorbed power law plus neutron 
star H atmosphere. The lower panel shows the ratio between the 
observed flux and the predicted model flux for each energy bin. See 
text and Table 1 for best-fit parameters. 



4.2. Orbital Phase-resolved Spectrum 

The ample number of photons (3280 within 2" of the 
pulsar) makes it possible to consider the continuum emis- 
sion from J1023+0038 from two portions of the orbit: 
(j>b = 0.0 — 0.5 (at the X-ray minimum) and 4>b = 0.5 — 1.0 
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Fig. 4. — Chandra ACIS-S3 1.5' X 1.5' image in J2000 coordinates 
centered on the X-ray counterpart to PSR J1023+0038 in the 0.3-8 
keV band. The image is binned at the intrinsic detector resolution 
of 0.5". The grey scale corresponds to counts increasing logarith- 
mically from white to black. The red arrow shows the direction 
of the pulsar proper motion (Deller et al. in preparation). Pixel 
randomization has been removed from the pipeline processing. 

the X-ray minimum the emission from the cooler por- 
tions of the polar cap(s) becomes more prominent. Al- 
ternatively, the differences in inferred temperature could 
be caused by the presence of both polar cap and "coro- 
nal" thermal plasma emission, especially if the coronal 
emission is also modulated at the binary period. De- 
tailed spin and orbital phase-resolved spectroscopic ob- 
servations are required to unambiguously disentangle the 
various spectral components. 

5. IMAGING ANALYSIS 

As noted in §1, evidence of past pulsar wind activa- 
tion and accretion episodes in PSR J1023+0338 may, 
in principle, be manifested through irregularities or dis- 
continuities in any nebular X-ray emission associated 
with the pulsar. The expected stand-off angular distance 
for a bow shock around a pulsar moving supersonically 
through the ambient medium can be estimated using the 
expre ssion 9 — n~ 1 / 2 /i~ 1 E 1 / 2 D~ 2 (Kargaltsc v k, Pavlovl 
120081 and references therein), where n is the mean den- 
sity of the ambient gas, fj, is the pulsar proper motion, 
E the pulsar spin-down luminosi ty, and D the distance . 
For a distance of D = 1.3 kpc (j Archibald et al.ll2009D . 

= 18 mas yr -1 , E — 5 x 10 34 ergs s _1 , and n k, 0.1 
cm~ 3 expected for warm n eutral ambient gas (see, e.g., 
IGhatteriee & Cordesl l2002h . we obtain a stand-off dis- 
tance of ~4". 

Figure 4 shows the image of the ACIS-S exposure cen- 
tered on PSR J1023+0038 in the 0.3 - 8 keV band. 
The binary is the brightest source in the image at the 
wavdetect-reported position of a = 10 h 23 m 47?68, S = 
+00°38'40"99. This differs from the optical position 
a = 10 h 23 m 4 7^67, 5 = +00°38'41'/ 2 given in the NO- 
MAD catalog (jZacharias et al.ll2004D by just +0.15" and 
—0.21" in right ascension and declination, respectively, 
well within the typical uncertainty in the absolute as- 




FlG. 5.— (Top) Chandra ACIS-S 0.3-8 keV image of the X-ray 
counterpart to PSR J1023+0038 from ObsID 11075. The color 
scale corresponds to counts increasing logarithmically from black 
to white. The green circles in all three images are of radius 2" 
and are centered on the X-ray source position. (Middle) Average 
of 100 simulated point spread functions of PSR J1023+0038 with 
ChaRT and MARX assuming the observing parameters of ObsID 
11075 and the spectrum of the pulsar as inferred from the formal 
spectral fits. The color scale is the same as in the top panel. (Bot- 
tom) The resulting difference image (observed minus simulated). 
The artifacts within l" are most likely caused by imperfections in 
the simulated PSF core. Beyond the PSF core, the residuals are 
consistent with being solely due to statistical fluctuations. 

trometry of ChandrdF\ It is evident that beyond ^2" of 
the pulsar, there is no significant excess of photons sur- 
rounding the pulsar, as might be expected from a pulsar 
wind nebula, nor in the direction opposite the pulsar 
proper motion as expected for a fast-moving pulsar. 

PSR B1957+20 {D = 2.5 kpc, E = 10 35 ergs s" 1 ) 
exhibits an X-ray tail of length ~16" in the direc- 
tion opposite the pulsar proper motion (Stapper s et al.l 
2003). Extracting the counts from a rectangular re- 
gion 20" x 5" trailing J1023+0038 yields no apprecia- 
ble excess of counts above the background count rate 
near the pulsar position, which is 2 x 10~ 6 counts s _1 

10 See http://cxc.harvard.edu/cal/ASPECT/celmon/ 
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Radial profile of the observed Chandra ACIS-S point 
spread function of PSR J1023+0038 (histogram) and the average 
of 100 simulated PSFs with ChaRT and MARX (red dashed line). 



arcsec -2 (0.3-8 keV). Assuming a representative power- 
law spectrum from the cur rent sample of X-ray PWNe 
(jKargaltsev fc Pavlovll2008h with r = 1.5, the conserva- 
tive limit on the wind nebula luminosity is ^3.6 x 10 29 
ergs s _1 for a distance of 1.3 kpc. 

To investigate the presence of small-scale (< 2") dif- 
fuse emission that could potentially arise due to a very 
compact pulsar wind nebula, we simulated a set of 100 
observations of PSR J1023+0038 with the ChaRT___ and 
MARX 4.5___ software packages, based on the parameters 
of the ACIS-S observation and the spectral shape of the 
X-ray emission as determined from the spectral fits (see 
§4). We subtracted the resulting average synthetic point 
spread function from the observed image and examined 
the significance of any residual emission relative to the 
background level (Figure 5). Within ~1" of the cen- 
ter, the difference image exhibits significant azimuthally 
asymmetric residuals (both positive and negative). One 
possible cause could be a slight offset between the loca- 
tion of the source centroid reported by wavdetect and 
the true source position. To test this possibility, we 
have simulated images with a range of displacements 
along the azimuthal direction in which the residuals are 
most pronounced in an effort to reduce the observed ar- 
tifacts. However, in all instances there is no improve- 
ment in the difference image. Therefore, it appears more 
likely that the artifacts are due to imperfections in the 
modeling of the High- Re solution Mirror Assembly optics 
(| Juda fc Karov ska 20100. Aside from the discrepancies 
in the central pixels of the PSF, the spatial distribution 
of events beyond ~1" of the pulsar position is consistent 
with that of a point source (Figure 6). 



11 The Chandra Ray Tracer, available a t 

http: / /exc .harvard. edu/sof t/ChaRT/cgi-bin/www-saosac . cgi 
Lz Available at http://space.mit.edu/cxc/marx/index.html 



6. A GEOMETRIC MODEL OF THE 
INTRABINARY SHOCK 

The morphology and spectral properties of the folded 
X-ray lightcurve in Figure 2 allow us to make impor- 
tant qualitative statements regarding the properties of 
the source of the non-thermal X-ray emission in the 
J1023+0038 system. For instance, the lack of significant 
accompanying spectral variability indicates that the X- 
ray variability is not due to photoelectric absorption. An 
enhanced 2Vh during eclipse phases would diminish the 
softest portion of the spectrum causing an apparent hard- 
ening of the emission. Instead, a uniform decline at all 
energies is observed. Thus, a more natural interpretation 
is obstruction of our view of the X-ray-emitting region by 
the companion star. Although it is difficult to pinpoint 
the exact start and end, the approximate duration of the 
X-ray eclipse (~0. 1-0.4 in orbital phase) implies that the 
source of the X-ray emission is much closer to the sec- 
ondary star (presumably near or at the inner Lagrangian 
[LI] point) than to the pulsar. For reference, for an as- 
sumed pulsar mass of 1.4 M Q , which implies an orbital 
inclination of 46°, the LI point is occulted for 32% of 
the orbit (</> 6 = 0.09 - 0.41). The rounded shape of the 
X-ray minimum and the gradual ingress and egress in- 
dicate that the region is gradually occulted, implying it 
is extended (of order the size of the star). The factor of 
~2-3 decline in the X-ray flux implies that at </>b ~ 0.25 
the secondary star blocks our line of sight to at least half 
of the emission region. 

To verify these statements, we have constructed a 
three-dimensional geometric model of the binary assum- 
ing a 1.4 M pulsar and i = 46° plus the orbital parame - 
ters measured from radio timing ([Archibald et alJl2009f ). 
Using this model, we generate lightcurves to compare 
against the Chandra data. For the entire range of orbital 
phases, we determined the fraction of the emission region 
under consideration that was occulted by the secondary 
star as viewed by the observer from a 46° angle. Based 
on this, we determined the observed flux relative to when 
the region is in full view. The resulting lightcurves were 
subsequently fitted to the data, with a flux scaling and 
DC offset as the only free parameters. 

One possibility is that the pulsar wind is interacting 
with material near the surface of the secondary star. If 
this is indeed the case, the shape of the X-ray-emitting 
region would roughly correspond to the shape of the side 
of the secondary star facing the pulsar (Figure 7). For 
this scenario, as illustrative examples, we consider both 
the case of a Roche-lobe filling star and one that is 90% 
the size of its Roche lobe. It is apparent from Figure 
8 that in both cases the resulting lightcurves (blue and 
magenta lines, respectively) approximate the shape of 
the observed modulations fairly well. 

An alternative origin of the observed X-rays is inter- 
action of the pulsar wind with matter issuing from the 
Roche- lobe fillin g companion through the LI point (see, 
e.g., Figure 2 in Bogd anov et al.ll2005h . As an approxi- 
mation to this outflow, we take a cylindrical volume with 
one of the circular bases centered on the LI point (Figure 
7). In this configuration, it is possible to account for an 
eclipse minimum occuring before <pb — 0.25 (as suggested 
hv trip data alheit at marginal significance) by consid- 



* See also http: //he™ -harvard ed^- j^m^s/^2010/HE^^^|^ht^ t j ff t from th direction of 

and http: //exc .harvard.edu/cal/Hrc/PSF/acis_psf _2010oct .html o o j 
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Fig. 7. — {Top) Schematic diagram illustrating the emission re- 
gions considered in the model. The orbital plane is in the plane of 
the page with the secondary star rotating counter-clockwise about 
the center of mass of the system. The labels mark the general 
directions to the observer at both conjunctions and quadratures. 
The blue area corresponds to X-ray emission from the face of a 
Roche lobe-filling secondary star. The magenta is for a compan- 
ion that is 90% the size of its Roche lobe. The green and red 
areas are for isotropically emitting cylindrical regions at inner La- 
grangian point (LI). Note the offset of the regions with respect to 
the semi-major axis in the direction opposite the sense of rotation 
of the star (dotted line). (Bottom) Same as above but with the 
orbital plane perpendicular to the plane of the page. The dashed 
line shows the line of sight from the LI point to the observer at 
superior conjunction (<f>i, = 0.25) assuming i = 46°. 

motion of the secondary star in the binary. Two rep- 
resentative examples of such a "stream" geometry that 
reproduce the rough shape of the Chandra lightcurve are 
shown in Figure 8 (green and red lines). 

We note that it is also very possible that both emis- 
sion from the stellar surface and from a stream at LI 
contribute to the observed non-thermal X-ray flux. How- 
ever, due to the limited photon statistics, substantially 
deeper observations are required to unambiguously de- 
termine the shock geometry by way of lightcurve model- 
ing. Nonetheless, Figure 8 illustrates that an intrabinary 
shock localized primarily near LI and/or at the face of 
the companion is a very plausible interpretation of the 
observed X-ray flux modulations. 

7. CONSTRAINING THE PHYSICS OF THE 
INTRABINARY SHOCK 

The purely geometric model employed above pays no 
regard to the actual physical mechanism responsible for 
the production of the observed X-ray emission. Never- 
theless, the geometric information obtained can be used 
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Fig. 8. — Resulting model lightcurves fitted to the observed 
Chandra lightcurve in the 0.3-8 keV band. The colors correspond 
to the shaded regions in Figure 7. The dashed lines mark the or- 
bital phase range at which the LI point is eclipsed by the secondary 
star for the Roche lobe filling case. In all instances, the nominal or- 
bital parameters of the J1023+0038 system are assumed. The bin 
with the maximum count rate at 4>\y = 0.5 — 0.6 has been excluded 
from the lightcurves fits. 



to place constraints on the physics of the shock. Given 
the small separation between the pulsar and the compan- 
ion (1.2 x 10 11 cm), the shock must occur in a relatively 
strong magnetic field. Therefore, synchrotron emission 
from accelerated charged particles is the likely X-ray pro- 
duction mechanism. The X-ray luminosity in the shock is 
dependent on both the postshock magnetic field strength 
and the wind magnetization factor a, defined as the ratio 
between the Poynting flux and particle flux (or, equiv- 
alently, the ratio of magnetic and particle energy densi- 
ties). Since a at the shock front is not known a priori, we 
must consider both a wind dominated by kinetic energy 
(a = 0.003 as seen in the Crab nebula) and one that is 
magnetically dom inated (a ^> 1). Follow ing the formal- 
ism described by lArons &: Tavanil (|1993f) , the magnetic 
field strength immediately upstream of the shock is given 

byBi = [a/ {l + a)]^ 2 (E /cf p r 2 ) 1/2 ■ Here f p is a geomet- 
ric factor that defines the fraction of the sky into which 
the pulsar wind is emitted, and r is the distance between 
the MSP and the shock. We take r ~ 1 x 10 11 cm, the 
distance from the MSP to Li assuming Mmsp = 1.4 M© 
and i — 46°. For an isotropic pulsar wind ( f p = 1) and 
E = 5 x 10 34 ergs s" 1 we obtain B x w 0.75 G (a = 0.003) 
and B\ « 13 G (er 3> 1). The implied magnetic field 
strength beyond the shock is then B 2 = 3-Bi ~ 2.1 G or 
B 2 ~ 39 G. Production of the observed e = 0.3 — 8 ~ 1 
keV photons via synchrotron emission requires relativis- 
tic electrons/positrons with Lorentz factors 7 = 2.4 x 
10 5 (e/ B2) 1 / 2 ~ 10 5 . The corresponding synchrotron ra- 
di ative loss time is t svnc h =5 .1 x lO^Bf) -1 - 10 - 700 
s (Rv bicki fc Lightmanlll979h . The geometric model pre- 



sented earlier indicates that the shock region is of order 
the diameter of the secondary star (r « 5 x 10 10 cm). 
From here, the residence times of the radiating parti- 
cles in the shock region are f now = c/3r 5.1 s (for 
cr = 0.003) and t flow = c/r « 1.7 s (for a > 1). 
The expected synchrotron luminosity from the shock 



can be estimated from / s h oc kAeL £ 



/synch /7/g 



,E, 



where / syn ch is the radiative efficiency of the synchrotron 
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emission, / 7 is the fraction of the intercepted spin-down 
energy flux that goes into accelerating electrons with 
Lorentz factor 7, / ge0 m is the fraction of the pulsar 
wind that interacts with matter from the companion, 
/shock is the portion of the total observed X-ray lumi- 
nosity that originates from the intrabinary shock, and 
L e is the synchrotron luminosity over the energy band 
Ae. For simplicity, we will assume Ae « 1 keV cen- 
tered at 1 keV. Since / syn ch = (1 + i S ynehAflow)~\ we 
obtain / syn ch ~ 0.007 and / syn ch ~ 0.16 for a = 0.003 
and a 3> 1, respectively. The companion star intercepts 
jgeom ~ 0.01 of the pulsar's outflow for the case of an 
isotropic wind. It is important to note, however, that the 
MSP outflow most likely exhibits significant anisotropy, 
with the bulk of t he wind emitted equatorially (as seen in 
the Crab pulsar: iHester et al.lfl995t lMichellll994[ ). Since 
the spin axis of the pulsar has likely aligned with the or- 
bital angula r momentum vector during the LMXB a ccre- 
tion phase ()Bhattacharva fe van den Heuvell IT991D . the 
wind of PSR J1023+0038 should thus be preferentially 
emitted in the orbital plane. In the extreme case of a 
sheet-like wind confined entirely to the orbital plane, the 
secondary star would receive / gcom w 0.07 of the pulsar 
wind's power. Based on this, we set 0.01 < f geom < 0.07. 

Using the count rate maxima predicted by the model 
lightcurves shown in Figure 8, the intrinsic synchrotron 
luminosity of the shock in the absence of eclipses is es- 
timated to be L £ 3 x 10 31 ergs s _1 for particles with 
7 w 10 5 in a band Ae = 1 kev centered on 1 keV. With 
/shock ~ 0.87 as deduced from the spectral fits, we ob- 
tain 0.93 < / 7 < 5.6 (a = 0.003) and 0.04 < / ? < 0.25 
(er ^> 1). For a = 0.003, the derived f 1 for an isotropic 
wind exceeds unity and even in the highly anisotropic 
case is unrealistically high. The same applies if we con- 
sider a factor of 2 — 3 larger emission region and stronger 
upstream magnetic field. Since it is likely that L < 1 
(e.g., / 7 = 0.04 for the Crab pulsar and nebula), the 
unphysical values of / 7 for a = 0.003 suggest that at 
the distance of the shock, the wind of PSR J1023+0038 
is magnetically dominated (a 3> 1). Furthermore, for 
u > 1, a value of f 1 comparable to that of the Crab 
pulsar requires significant anisotropy in the wind. 

It is interesting to highlight the apparent enhancement 
of emission at eclipse egress (</>& = 0.5 — 0.6) since it likely 
arises due to D oppler beaming in the direction of the 
post-shock flow (|Arons fe Tavanilll993l ). Similar, albeit 
much less statistically significant features have been seen 
in other eclipsing MSP syste ms (see, e.g., IStappers et all 
12001 IBogdanov et"all 120051) . For PSR J1023+0038, the 
implied flow direction is opposite to the sense of motion 
of the secondary star in its orbit, indicating that the 
combined effects of the gravitational, pulsar wind pres- 
sure, and Coriolis forces together with the outflow veloc- 
ity relative to the orbital velocity, conspire to force the 
outflowing gas towards the trailing e dge of the eclipse 
(see, e.g.. iTavani fe Brookshawlll99lD . This effect can 
also account for the absence of a matching X-ray flux 
increase at ingress as it results in a weaker outflow to- 
wards the leading edge of the eclipse. Within the avail- 
able photon statistics, the factor of 1.2 — 1.3 increase 
in brightness at cj>b = 0.5 — 0.6 relative to the flux out- 
side of eclipse (fa ~ 0.6 — 0.9) is in agreement w ith the 
range 1.3 — 2.2 predicted bv lArons fe Tavanl (|1993f ) based 
on the expected range of post-shock wind flow velocities 



(c/3-( C /3)V3). 

8. CONCLUSION 

We have presented Chandra ACIS-S observations of 
the PSR J1023+0038 binary system. The data defini- 
tively confirm the X-ray flux variability at the binary or- 
bital period. A simple geometrical model can account for 
the general shape of the X-ray lightcurve of the binary, 
indicating that the observed variability is primarily due 
to a partial geometric eclipse of the intrabinary shock by 
the secondary star. The model suggests that the shock is 
likely localized near LI and/or is formed by the interac- 
tion of the pulsar wind and matter at the surface of the 
secondary star. The energetics of the shock emission fa- 
vor a magnetically dominated (a> 1) pulsar wind at the 
shock front. This is not surprising, given that the shock 
occurs in a ^10 4 times stronger magnetic field than in 
the case of the Crab p ulsar (where the wind is d ominated 
by kinetic energy; see lKennel fe CoronitilH984T ). The X- 
ray data also favor a significantly anisotropic pulsar wind 
that is preferentially emitted in the orbital plane, which, 
in turn, requires close alignment of the pulsar spin and 
orbital angular momentum expected from binary 

evolution. 

Com bined with the study presented in I Archibald et all 
( 2010), the spectral continuum suggests a dominant non- 
thermal component and the presence of a fainter ther- 
mal component that, at least in part, originates from 
the hot polar caps of the pulsar. In this sense, in terms 
of its X-ray properties, the pulsar itself appears to be 
similar to the maj ority of "recycled" MSPs detected 
in X-rays to date (iZavlinl 120061: IBogdanov et al.| [ 20061 : 
IBogdanov fe Grindlavl 120091 : IBogdanov et al.l 120111 1. A 
portion of the thermal X-ray emission may also come 
from an optically-thin thermal plasma, possibly the same 
plasma responsible for the radio eclipses. 

We find no evidence for a wind nebula associated with 
the pulsar. The shape, luminosity, and spectrum of a 
PWN depend on the angular distribution, magnetization 
and energy spectrum of the unshocked wind, as well as 
on the pulsar velocity and the properties of the ambient 
medium. Therefore, the absence of a nebula down to a 
fairly deep limit of <3.6 x 10 29 ergs s _1 , corresponding to 
<7 x 10~ 6 of the pulsar's E, implies that the combination 
of low ambient density, moderate space velocity, and/or 
unfavorable wind geometry are insufficient to drive an 
X-ray-luminous bow shock. 

Owing to its proximity, PSR J1023+0038 is one of the 
best-suited targets for studies of these rare "transition" 
systems and thus warrants further X-ray observations, 
especially in combination w ith optical spectr oscopy and 
7-ray variability study (see lTam et alJl2010L for the re- 
cent detection of this system with the Fermi Large Area 
Telescope). Constraining the properties of the intrabi- 
nary shock in J1023+0038 and similar systems has im- 
portant implications for understanding the X-ray prop- 
erties of accreting MSPs in qu iescence, such as the pro- 
totyp ical SAX J1808.4-3658 (|Wijnands fe van der Klisl 
1998). It has been suggested that the non-thermal X-ray 
radiation observed from certain quiescent low-mass X-ray 
binaries is also due to an active rotation-powered wind 
but this cannot be d irectly confirmed due to the absence 
of radio pulsations (|Campana et al.ll2004D . Alternative 
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scenarios such as o n-going, low-level accretion (see, e.g., 
iGarcia et al.ll2001 . and references therein) or a propeller 
regime (|Wiinandsl2003l and references therein) have also 
been proposed as explanations for the non-thermal emis- 
sion. An X-ray spectroscopic and variability study anal- 
ogous to that presented above could reveal the signature 
of an interacting rotation-powered pulsar wind in such 
systems, thereby providing further evidence for the evo- 
lutionary relation of LMXBs and MSPs. 
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